Muon spin rotation study of the intercalated graphite 
superconductor CaCg at low temperatures 

F. Muranyi\ M. Bendele^•^ R. Khasanov^, Z. 
Guguchia^, A. Shengelaya^, C. Baines^ and H. Keller^ 

^ Physik-Institut der Universitdt Zurich, 
Winterthurerstrasse 190, CH-8057 Zurich, Switzerland 
"^Laboratory for Muon Spin Spectroscopy, 
Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland 
^Department of Physics, Tbilisi State University, 
Chavchavadze 3, GE-0128 Tbilisi, Georgia 
(Dated: August 29, 2011) 

Abstract 

Muon spin rotation (^uSR) experiments were performed on the intercalated graphite CaCe in the 
normal and superconducting state down to 20 mK. In addition, AC magnetization measurements 
were carried out resulting in an anisotropic upper critical field ifc2) from which the coherence 
lengths i^gft (0) = 36.3(1.5) nm and (0) = 4.3(7) nm were estimated. The anisotropy parameter 
■^^ = ff^|/if^2 increases monotonically with decreasing temperature. A single isotropic s-wave 
description of superconductivity cannot account for this behaviour. From magnetic field dependent 
/iSR experiments the absolute value of the in-plane magnetic penetretion depth \ab{^) = 78(3) nm 
was determined. The temperature dependence of the superfiuid density (T) is slightly better 
described by a two-gap than a single-gap model. 
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I. INTRODUCTION 



The field of grapliite intercalation compounds (GICs) gained attention after the discov- 
ery of the superconductor CaCe with a rather high value of the superconducting transition 
temperature ~ 11.5 K [1]. Superconductivity in GICs was first reported in the potassium- 
graphite compound KCg with ^ 0.14 K |2|. Until the discovery of the superconductor 



CaCe the hig. 
ent pressure 



lest Tc was observed in KTli 5C5 with ~ 2.7 K, synthesized under ambi- 
3|. However, high pressure synthesis was found to increase Tc up to 5 K in 



metastable compounds such as NaC2 and KC3 . According to experimental and theo- 
retical work CaCg can be described as a classical BCS superconductor with a single isotropic 

7|. The upper critical field ifc2 shows a remarkable anisotropy with 



gap Aq 1.7 meV 



all 



zero temperature values /io-f^c2(0) — 0-3 T and /io-f^c2'(0) ~ 1.9 T, for the external field H 



parallel to the c-axis or parallel to the afe-plane, respectively jsl. A recent ARPES study 
indicated the existence of a possible second superconducting gap jo] with a small zero tem- 



perature value Ao,2 — 0.2(2) meV. Tunneling experiments [lO!] gave strong indications for the 
existense of a superconducting anisotropic s-wave gap in CaCg, providing different zero tem- 
perature gap values for changing current injection directions, Aq ^ 1.7 meV and Aq'' — 1.44 
meV. The muon-spin rotation (yuSR) technique is a powerful method to characterise the su- 



perconducting gap symmetry in superconductors 



ll| . /iSR measurements down to very low 



temperatures may allow to access a possible second small gap, which should manifest itself 
in the low temperature superfluid density in terms of an inflection point [llj. A recent /iSR 
study [l2| supports a single gap isotropic s-wave description of superconductivity in CaCe, 
although the low temperature region was not studied. In this work we present extended 
measurements down to 20 mK and investigate possible order parameter symmetries (single- 
gap isotropic s-wave, single-gap anisotropic s-wave, two-gap isotropic s-wave) by means of 
liSR experiments. 

II. EXPERIMENTAL RESULTS 

CaCe samples were prepared from Highly Oriented Pyrolytic Graphite (HOPG, Structure 
Probe Inc., SPI-2 Grade) and calcium metal (Sigma Aldrich, 99.99% purity) by means of 



the molten alloy method. The detailed preparation procedure is described elsewhere 13 |. 
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In order to reduce the melting point of the calcium, lithium (Sigma Aldrich, 99.9% purity) 
was added. The atomic ratio 3:1 of calcium and lithium was measured in a stainless steel 
reactor under argon atmosphere. The molten alloy was cleaned in the glove-box. After 
immersing the graphite pieces the reactor was tightly sealed. Intercalation at 350 °C was 
carried out in an industrial furnace for 30 days. The resulting sample's color is gold and 
the surface possesses a shiny metallic lustre. The final sample is a superconductor with an 
onset temperature — 11.2 K, in reasonable agreement with earlier studies [l, 6, 7|. 

To characterise our sample we carried out detailed AC magnetization measurements pre- 
sented in Fig. 1. The AC component of the applied magnetic field was 1 mT 
at 1 kHz. Figure la shows typical AC magnetization curves at 80 mT (real M'{T) and 
imaginary M"{T) part), the straight lines indicate the onset criterion to determine T^. The 
temperature dependence of the upper critical field components (-f^lk) and {H\\ab) 
are displayed in Fig. lb. The values of the and were determined from AC magneti- 
zation using the procedure illustrated in Fig. la. The upper critical field is anisotropic and 

n n 

follows the temperature dependence as reported in Refs. [8[ and p^. In the framework of 
the anisotropic Ginzburg-Landau theory the upper critical field for H\\c and H\\ab is given 
by 



A^o^ct (0) ^ 



2vrC. (0) C (0) ' 

/^o^c2(0) = -J^, (1) 

where $o = h/2e = 2.07 ■ 10^^^ Tm^ is the flux quantum, and ^^^^ are the corresponding 
coherence length components at zero temperature. With the linearly extrapolated values of 
the upper critical field at zero temperature, fiQH^2 i.^) = 2.1(3) T and yUQif^2 (0) = 0.25(2) 
T, and the equation above we obtain (0) = 36.3(1.5) nm and C,^ (0) = 4.3(7) nm. These 
values are comparable to earlier reports such as resistivity measurements (^^^ (0) = 29.0 nm 
md (0) = 5.7 nm) j8|, or susceptibility studies (■Cat (0) = 35 nm and (0) = 13 nm) 
15| . Remarkably, the upper critical field H^2 shows a positive curvature in the studied 
temperature region which was not observed in previous investigations. A similar positive 
curvature was observed near in MgB2 [16] where it was explained by a two-gap model. 
However, we should keep in mind that other explanations are also possible. Complementary 
experiments, such as studies of the superfiuid density, as performed in this work, are nec- 
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essary to draw definite conclusions. The temperature dependence of tlie upper critical field 
anisotropy 7^^^ = H'^2/H^2 is shown in the inset of Fig. lb. We interpolated the measured 
upper critical field values with a third order polynomial and determined the ratio '-jjj. The 
anisotropy paremeter '-yjj increases monotonically with decreasing temperature, showing a 
similar temperature dependence as with MgB2 17|]. Note that a single gap isotropic s-wave 
description of superconductivity cannot account for this behaviour of 7^:^ (T). 

/xSR experiments were carried out at the vrMS beamline using the GPS (General Purpose 
Spectrometer) and the LTF (Low Temperature Facility) spectrometers at the Paul Scherrer 
Institute PSI Villigen, Switzerland. Six pieces of CaCg forming an area of ~ 10x14 mm^ 
were used in the experiment. In the intercalation process the Ca atoms penetrate from the 
side of the graphite sample and diffuse along the a6-plane. Reducing the sample size in the 
a6-direction favors Ca diffusion, and consequently leads to a higher sample quality than just 
using one large piece. For the LTF experiments the samples were glued on a silver plate with 
the help of Apiezon N grease and covered by silver coated polyester foil. In contrast, for the 
GPS experiments the samples were supported between the two arms of a specially designed 
fork by means of kapton tape. This technique allows only the incoming muons that stop in 
the sample to be counted, and therefore the background signal is reduced. The magnetic 
field (80 mT) was applied along the crystallographic c-direction at 15 K well above T^. Then 
the sample was cooled down to the lowest temperature (20 mK for the LTF spectrometer), 
and /iSR spectra were recorded with increasing temperature. The same procedure was 
applied in the GPS spectrometer, down to the lowest temperature of 1.6 K available in this 
device. The initial spin polarization of the implanted muons is perpendicular to the applied 
magnetic field H. The /iSR experiments were performed at H > H^^ so the sample was 
in the vortex state. Magnetic field dependent measurements were performed in the GPS 
spectrometer. At each magnetic field H the sample was warmed up to 15 K (well above Tc), 
a /iSR time spectrum was recorded, then it was cooled down to 1.6 K, and another /iSR 
time spectrum was taken. Before and after the /iSR experiments the same piece of 
CaCe was measured in the AC susceptometer, comparing the two measurements 
we observed no difference. 

In the vortex state of a type II superconductor the muons probe the local magnetic 
field distribution P{B) due to the vortex lattice [l^. From the second moment (Ai?^) 
of P{B) the in-plane magnetic penetration depth \ab can be extracted according to the 
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relation (AS^) = C\~^ {C is a field dependent quantity) from which the superfiuid density 
Ps oc can be determined Thus, pSR provides a direct method to measure the 
superfiuid density in the mixed state of a type II superconductor. Figure 2a shows the /iSR 
asymmetries A{t) recorded in LTF in zero magnetic field above and at 20 mK. The two 
time spectra overlap, indicating the absence of any kind of magnetic order in the sample and 
in the sample holder. The time dependence of the asymmetry can be described by the static 
Kubo-Toyabe formula, showing only the presence of randomly distributed nuclear magnetic 
moments, but no realization of magnetism. Figure 2b shows the /iSR asymmetries taken at 
15 K well above Tc ^ 11.5 K and at 20 mK in a field of 80 mT. A clear damping of the pSR 
signal at 20 mK due to the presence of the vortex lattice is visible. For a detailed description 
of /iSR studies of the vortex state in type II superconductors see, e.g. [isl. The pSR time 
spectra were fitted to the following expression: 

A^^it) = Asce('^^«*)'/2cos(7^i?sc + 0) 

+A^,eMy^cos{^^B^, + cP) (2) 

+AAg COs(7^5Ag + 0) 



Here the indices SC, bg and Ag denote the sample (superconductor), the background 
arising from non superconducting parts of the sample and the non-relaxing silver back- 
ground, respectively. A denotes the initial asymmetry, a is the Gaussian relaxation rate, 
7^ = 27r ■ 135.5342 MHz/T is the muon gyromagnetic ratio, B is the internal magnetic field, 
and is the initial phase of the moun-spin ensemble. A set of /iSR data were fitted simul- 
taneously with y4sc, ^bg, ^Ag, -Bbg, -BAg, (Xbg, and (p as common parameters, and Bsc and 
(Tsc as free parameters for each temperature. The same expression was used to analyse the 
data recorded with the GPS spectrometer (where no silver background signal is present). 
To analyze the magnetic field dependent /iSR measurements -Bsc was a free parameter in 
Eq. (2) since it depends on the applied field. 

The analysis of the field dependence of the measured /iSR relaxation rate (Xsc oc {AB"^) 
allows to estimate the absolute value of the in-plane magnetic penetration depth Aab, as- 
suming that CaCe can be described as a single isotropic s-wave gap superconductor. The 
measured values of asc are plotted in Fig. 3. To analyze our data we used the formula 
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developed by Brandt |19| : 



$0 



CISC ^ 0.1727„-^ (I- — 



271 



B 



B. 



c2 



B 



1 + 1.21 I- ^ — 



B: 



c2 



A 



ab ■ 



(3) 



Although the formula was derived for high k, type II superconductors {k > 5), it was 
suggested to work also for low values of k at mag netic fields ^jBjBl^ > 0.5 jl2|. From 
the measured fiSR relaxation rates asc and Eq. (3) one obtains Xab (0) = 78(3) nm and 
5^2(0) = 170(5) mT. The solid line in Fig. 3 represents the corresponding fit (the point at 
200 mT was not included in the fit). To estimate the Ginzburg-Landau parameter k we use 
our values of ^^j, (0) = 36.3(1.5) nm and A^;, (0) = 78(3) nm giving k ^ 2.1. It is worth to 
mention that the value of B^^ is smaller than that estimated from the AC magnetization 
measurements (5^2(0) — 250 mT), since AC magnetization measurements resulted 
more data points to determine the upper critical field value we tctke 250 mT as 
5^2(0) for the further analysis. 

It is generally assumed that the superfluid density pg is related to the in-plane magnetic 
penetration depth Aab by the simple relation pg oc \~^. However, for magnetic fields close 
to Bc2 this proportionality must be corrected because the order parameter ip (r) is reduced 
due to the overlapping vortices. In this case the spatial average of the superfluid density is 
given by |l2|: 



p, (T) ^ (1^ (r) \')\-' (T) ^ (1 - B/B:,iT)) (T) , (4) 

where (...) means the spatial average. Values of were calculated using Eq. (3) and for 
the values of B^2 we used the results of our AC magnetization measurements (see Fig. lb). 
The temperature dependence of the superfluid density (T) determined at 80 mT is plotted 
in Fig. 4. We analyzed p^ (T) using three different models: i) single-gap isotropic s-wave ii) 
single-gap anisotropic s-wave, and in) two-gap isotropic s-wave. To calculate the tempera- 
ture dependence of the magnetic penetration depth within the local approximation (A ^ ,^) 
we applied the following formula 20|, l2l| : 



pAt) = Ps{o) 1 + 



TT 



2lT 



dl 
dE 



E dEdcp 



(5) 



where pg (0) is the zero temperature value of the superfluid density, / = 
[1 + exp {{E — Ep) / k-sT)]'^ is the Fermi function, is the angle along the Fermi surface, 



and A (T, 0) = AqH (T/Tc) g {(p). The temperature dependence of the gap is expressed by 
h (T/Tc) = tanh (1.82 [1.018 {TJT - 1)]°-^^) For the isotropic s-wave gap i) g"" {(f)) = 1 
and for the anisotropic s-wave gap ii) g^^'^ (0) = (1 + a cos 40) / (1 + a), where a denotes 
the anisotropy of the gap. The two- gap fit was calculated using the so called a-model and 
assuming that the total superfiuid density is the sum of the two components 



21 



22 : 



Ps iT) = Ps . ^ + (1 - w) , , . (6) 

V P(0,Ao,i) p(0,Ao,2)/ 

Here Aq^i and Ao,2 are the zero temperature values of the larger and the smaller gap, 
and u; (0 < w < 1) is a weighting factor representing the relative contribution of the larger 
gap to Pg. The results obtained for the different models are plotted in Fig. 4. The black 
line represents the result for the single s-wave gap model i) with Pg (0) = 96.5(5) pm~^, 
Ao = 0.79(1) meV, = 6.80(2) K, and 2Ao/A;bTc = 2.70(5), describing CaCg as weak 
coupling BCS superconductor. The reduction of agrees well with our AC magnetization 
measurements for 80 mT. The blue, dotted line represents the anisotropic s-wave gap analysis 
ii) with Aq = 0.81(1) meV and a ^ 0.37 for the gap anisotropy. These values do not agree 
with the reported ones from tunneling experiments (Aq ^1.7 meV and Aq* ~ 1.44 meV) 



10|. The two-gap analysis Hi) (red, dashed line) yields Ao,i = 0.85(2) meV for the larger 



gap and Ao,2 = 0.23(1) meV for the smaller one. Only ~ 9.4(3) % of the superfiuid density 
is associated with the smaller gap. The smaller gap agrees with Ao,2 = 0.2(2) meV reported 
in js]. Our larger gap value is significantly smaller (Ao,i = 1.9(2) meV). It is worth to 
mention that the applicability of the two gap model at this magnetic field (80 
mT) is questionable. The gap magnitude scales with the upper critical field, 
(zTl) ~ i^' resulting i?^2 ~ 12 mT which is much smaller than the applied 
magnetic field. Although the best agreement is found for the two-gap scenario, no final 
conclusion can be drawn from this analysis. All three scenarios describe well the observec 
temperature dependence of the superfiuid density. For a direct comparison with Ref. [12 1 
we omitted the values of p^ (T) below 2 K. In this case the single s-wave gap fit yields 
Pg(O) = 91.6(9) pm-2, Ao = 0.84(1) meV, = 6.75(2) K, and 2Ao/kBT^ = 2.90(5). 
The values presented in Ref. [12|] are: (0) = 85.1(3) pm~^, Aq = 0.868(5) meV, and 
2Ao//cbT'c = 3.6(1). The higher applied magnetic field (120 mT in Ref. p^]) may account 
for the differences. Here, one should point out that our low temperature data are essential 
to determine the zero temperature values of Pg and A. 
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III. CONCLUSIONS 



In conlusion, we carried out extensive AC magnetization measurements to map the tem- 
perature dependence of the upper critical field and CaCg and to evaluate the 



coherence length. The values of ^^f, (0) = 36^3 
agreement with those reported previously 



1.5) nm and (0) = 4.3(7) nm are in good 



15|. We found that the upper critical field 



anisotropy ^y^j = H^2/H^2 increases with decreasing temperature. From magnetic field de- 
pendent /iSR experiments the absolute value of the in-plane magnetic penetration depth 
was determined to be \ab — 78(3) nm, in agreement with previously reported values [7, 112I . 
Furthermore, low temperature /iSR experiments were performed in order to map the whole 
temperature dependence of the superfiuid density Pg (T). We analyzed the temperature 
dependence of Pg with three different models: i) single-gap isotropic s-wave, ii) single-gap 
anisotropic s-wave, and in) two-gap isotropic s-wave. All models describe the measured 
pSR data almost equally well, although a slightly better agreement was achieved using the 
two-gap model. 

This work was partly supported by the Swiss National Science Foundation (SCOPES 
grant No. IZ73Z0_128242). The pSR experiments were performed at the Swiss Muon Source, 
Paul Scherrer Institut, Villigen, Switzerland. 
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FIG. 1: a) AC magnetization curves of CaCe taken at fi^H = 80 niT and niT. Af and M" denote 
the real and the imaginary part of the AC magnetization, respectively. The arrow indicates the 
onset criterion to determine Tc{H) (crossing of the two hnes). b) Temperature dependence of the 
upper critical fields ^^'^ the external magnetic field H\\c and H\\ab, respectively, as 

determined by AC magnetization. The inset shows the calculated upper critical field anisotropy 
'jfj = -f^c2/^c2 ^ function of temperature. 




FIG. 2: (color online) Measured /iSR asymmetries A{t) in CaCg. a) Asymmetry in zero magnetic 
field (ZF). b) Asymmetry in 80 mT transverse field (TF). The black data points were taken at 15 
K (well above Tc) and the red ones at 20 mK. 
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FIG. 3: Field dependence of the Gaussian relaxation rate use in CaCg measured at T = 1.6 K. 
The solid line corresponds to the calculated curve using Eq. (3) (see text). 
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FIG. 4: (color online) Temperature dependence of the normalized superfluid density Ps(^)/Ps(0) 
CaCe measured at 80 mT. The inset shows the low temperature region. The black line corresponds 
to a single-gap isotropic s-wave, the blue, dotted line to a single-gap anisotropic s-wave, and the 
red, dashed line to a two-gap isotropic s-wave description of Ps(T). For a detailed explanation see 
text. 
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